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1.  INTRODUCTION 


In  many  industrial  processes  and  research  related  efforts,  the  rapid 
and  reliable  disposal  of  large  quantities  of  gaseous  hydrogen  is  of  con- 
siderable concern.  Currently,  two  techniques  are  in  general  use.  The 
Kennedy  Space  Center  in  Florida  employs  a series  of  small  separate  ignition 
devices  to  burn  hydrogen  in  an  area  called  the  "hydrogen  burn  pond."  This 
technique  was  used  to  dispose  of  excess  hydrogen  employed  in  the  Apollo 
program.  The  Rocketdyne  test  facility  at  Canoga  Park,  California,  employs 
rather  high,  large  diameter  vent  stacks  with  pilot  flames  to  dispose  of  the 
hydrogen  used  in  testing  space  shuttle  components.  An  alternate  approach 
for  the  rapid  disposal  of  hydrogen  also  appears  possible.  Large  quantities 
of  hydrogen  may  be  released  in  "bursts"  or  "bubbles"  which  are  allowed  to 
mix  with  air  and  then  are  ignited  at  some  predetermined  time  or  altitude. 

This  report  explores  the  various  problems  associated  with  the  release  and 
subsequent  Ignition  of  these  "bubbles"  of  hydrogen. 

The  effort  has  been  divided  into  three  distinct  phases.  The  first 
phase  Involves  the  formulating  of  the  equations  of  motion  and  an  analysis 
of  the  vertical  and  horizontal  movement  of  the  bubbles.  The  predicted 
values  are  compared  with  the  data  of  Gorev,  Gusev,  and  Troshln*  and  the 
GEST  experiments^.  Next,  the  ignition  process  is  investigated  and  compared 
to  the  data  of  Gorev  et  al.  for  small  hydrogen  bubbles.  The  agreement 
generally  is  quite  good.  Finally,  the  formulation  for  the  infrared  radiation 
estimates  in  the  2.8  y band  of  water  is  discussed  and  predictions  are  com- 
pared with  the  GEST  experiments.  The  technique  for  estimating  the  radiant 
intensity  is  employed  to  predict  the  magnitude  of  radiation  to  be  expected 
from  the  burning  "bubbles"  of  hydrogen. 


2.  SIMPLIFIED  THEORY  OF  RISING  THERMALS 
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The  motion  of  a buoyant  thermal  has  been  explored  in  the  simplified 
models  of  Wang^'^  and  Shul  and  Weyl^  and  can  be  applied  directly  to  a 
bubble  of  gas  released  in  a stable  atmosphere.  The  momentum  equation  for 
a spherical  gas  bubble  of  radius  r , volume  V , and  density  p less  than 
of  the  ambient  medium  is: 

+ - (p„-  p)gv  (1) 

where  U is  the  velocity  of  the  center  of  mass  and  - .Sp^^V  is  the 
virtual  mass.  Form  drag  is  neglected  since  it  is  small  compared  with  the 
effective  drag  due  to  entrainment.  The  mass  entrainment  coefficient  o is 
defined  by  the  equation  of  mass  conservation: 


The  change  in  energy  of  the  bubble  is  due  primarily  to  the  internal 
energy  of  entrained  air,  which  is  large  compared  to  the  change  in  kinetic 
energy,  pressure  volume  work,  or  the  potential  energy  due  to  gravity.  Con- 
servation of  internal  energies  together  with  Eq.  2 yields: 

g ■ au  (3) 

Following  the  procedure  outlined  in  Shul  and  Weyl’  and  using  the 
dimensionless  variables,  x ■ " U/[2rj/2o(l  - , and 

T ■ t/(2rg/aAg)*^^  , an  expression  is  obtained  which  relates  the  expansion 
ratio  x to  the  dimensionless  time  T after  release: 

T - - 'Y  + 7 juj  x"  - AAx  - (1  - 4A)J  ‘ (4) 

where  Ug  is  the  normalized  initial  velocity  and 

A = 2(p^  - Po)/3P»  (5) 
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The  height  of  the  rising  bubble  z is  directly  proportional  to  the 
change  in  radius  (see  Figure  1); 

z “ ^r-rjj/o  ■ rjCx-D/a  (6) 

The  velocity  for  a given  x can  be  determined  from  Eq.  3: 


a 


— (2rj/aAg) 


dr 


(7) 


where  dx/dT  can  be  found  from  Eq.  4. 

2.1  Effect  of  Crosswind 

The  effect  of  a crosswind  on  the  rise  of  a hydrogen  bubble  has  also 
been  examined.  Form  drag  is  again  neglected  and  the  momentum  equation 
becomes : 


v(47rr^p„aVu^)+ ^ 


(8) 


where  dM^/dt  is  determined  from  Eq.  2,  U is  the  upward  velocity,  and  v 
is  the  horizontal  velocity  of  the  bubble  relative  to  the  air.  The  small 
change  in  upward  velocity  U associated  with  increased  entrainment  in 
crosswinds  has  been  ignored. 

Eq.  8 Is  Integrated  to  obtain: 


where  v starts  with  a negative  value  and  goes  to  zero  as  the  bubble  approaches 
wlndspeed.  Eq.  9 can  be  integrated  numerically  until  v has  dropped  to  5Z  of 
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its  original  value,  after  which  it  is  assumed  that  the  crosswind  no  longer 
affects  the  rise  of  the  bubble.  A bubble  20  meters  in  diameter  reaches  a 
crosswind  velocity  of  10  mph  in  only  .5  second  compared  to  .2  second  for  a 
1 meter  diameter  bubble.  If  the  bubble  has  an  initial  upward  velocity,  it 
reaches  crosswind  velocity  in  a slightly  longer  time.  Since  entrainment 
due  to  crosswind  is  of  such  short  duration,  the  major  effect  of  a crosswind 
is  to  move  the  bubble  downstream  by  a distance  x according  to  the  relation: 


The  major  effects  of  a 10  mph  crosswind  are  shown  in  Figure  2 for  bubbles 
of  1 meter  and  20  meters  initial  diameter. 

2.2  Measured  Rise  Times 

Gorev,  Gusev,  and  Troshin*  have  measured  the  rise  times  of  hydrogen 
bubbles  with  initial  volumes  of  approximately  70  liters.  Their  data  are 
presented  in  Figure  3. 

It  should  be  noted  that  Gorev  et  al.  define  their  entrainment  coeffi- 
cient a'  in  terms  of  the  velocity  of  the  top  of  the  bubble  rather  than 
the  center.  The  relation  between  their  empirically  determined  o'  and  the 
a defined  by  Eq.  2 is  therefore: 

a - a'/(l  -a') 

The  two  curves  shown  in  Figure  3 are  calculated  from  the  model  discussed 
above  for  a .41  (a'  ■ .29)  and  dimensionless  initial  velocities  of 
Fj  * 0 and  F,  » 4 where  Fq  - Uq/  '\Jr^g  . Data  and  prediction  are  in 
good  agreement. 

The  GEST  experiments^’®  form  another  set  of  empirical  data  in  which  the 
rise  times  of  thermals  have  been  measured.  In  the  GEST  experiments,  a 
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Figure  3.  Rise  time  of  H2  bubbles  — from  Refer' 


10  meter  diameter  mylar  balloon  filled  with  a mixture  of  1 part  CH^  to 
2 parts  O2  was  ignited  at  an  altitude  of  about  ISO  ft.  The  height  from  the 
ground  to  the  center  of  the  rising  toroidal  fireball  was  measured  as  a 
function  of  time  after  ignition.  The  data  are  shown  In  Figure  A.  The 
theoretical  rise  rate  calculated  from  the  model  agrees  with  the  GEST  data 
if  the  entrainment  coefficient  Is  taken  to  be  as  .1-.12  Instead  of  .41 
used  by  Gorev  et  al.  Predictions  for  these  values  of  a are  shown  In 
Figure  4.  Infrared  radiation  time  histories  were  also  measured  In  the  GEST 
experiments.  It  will  be  shown  In  a later  section  that  a value  of 
a “ .1-.12  is  also  required  to  correlate  the  measured  and  predicted 
radiation  levels  as  a function  of  time.  Still  another  set  of  data  on  the 
rise  times  of  thermals  caused  by  nuclear  explosions  may  be  found  In 
Reference  5 (also  see  References  7 and  8).  Apparently,  an  entrainment 
coefficient  of  .25  Is  required  to  represent  these  very  large  thermals. 

2.3  Choice  of  Entrainment  Coefficient  for  Large  H2  Bubbles 

From  the  three  sets  of  data  In  the  previous  section  It  appears  that 
the  simplified  model  described  earlier  is  adequate  to  predict  the  rise  times 
of  buoyant  gas  bubbles  If  a suitable  entrainment  coefficient  Is  chosen. 
Assuming  that  the  observed  variations  In  a are  not  due  simply  to  experi- 
mental error.  It  appears  that  small  hydrogen  bubbles  ( ~ 70  liters)  entrain 
air  more  rapidly  than  either  very  la^ge  thermals  ( - 4 x 10^  cubic  meters) 
from  nuclear  blasts  or  large  gas  clouds  from  the  combustion  ^ducts  of 
methane-oxygen  bubbles  which  are  much  hotter  than  the  ambient  air. 

At  present  it  Is  not  obvious  which  choice  of  a Is  appropriate  for 
large  hydrogen  bubbles  with  volumes  from  .5  x 10^  to  4 x 10^  cubic  meters. 
If  one  assumes  that  the  initial  size  of  the  bubble  does  not  affect  the 
entrainment  rate,  then  larger  hydrogen  bubbles  must  behave  like  their  smaller 
counterparts  and  a will  be  close  to  .41.  Operating  under  this  assumption, 
an  a of  .41  was  used  In  all  the  hydrogen  bubble  calculations  presented  In 
this  paper.  It  Is  Important  to  note,  however,  that  after  ignition  the 
hydrogen  bubbles  may  behave  similarly  to  the  thermals,  requiring  a smaller 
a In  that  domain. 
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Rise  times  of  methane-oxygen  combustion  products. 


3.  MODELING  OF  THE  COMBUSTION  PROCESS 


In  addition  to  data  on  the  motion  of  small  hydrogen  bubbles,  Gorev  et 
al.  have  measured  the  ratio  ir  of  the  volume  of  the  combustion  products 
V2  to  the  bubble  volume  before  ignition,  V , as  a function  of  the  height 
at  which  ignition  occurs.  The  data  is  presented  in  Figure  5 where  the  solid 
line  is  an  empirical  expression  for  the  expansion  data  as  a function  of  the 
nondimenslonal  height  ( - h/r,  and  is  given  by 


x 


25  - '’♦.2C 

(.3C  + 1)* 


+ 1 


(10) 


If  it  is  assumed  that  the  expansion  represented  by  Eq.  10  is  independent 
of  initial  bubble  size,  an  expanded  volume  V2  ~ xV  can  be  assigned  to 
the  bubble  after  ignition  at  any  altitude  ( . Parameters  of  the  combustion 
process  such  as  percentage  of  hydrogen  burned,  degree  of  mixing,  and  amount 
of  air  entrained  after  ignition  can  be  obtained  from  models  which  maintain 
the  observed  expansion. 

I 

The  overall  chemical ^change  within  the  bubble  can  be  written  as 

1 H2  + y(.21  0-  + .79  Nj)  (w6)h20 

+ ^1  - wB^H2  + (.21y  - .5we)02  + (•79y)N2  (11) 


where  y is  the  overall  ratio  of  air  to  hydrogen  within  the  bubble  as 
determined  in  Section  2.0.  Due  to  incomplete  mixing  on  a molecular  scale, 
both  unreacted  oxygen  and  unreacted  hydrogen  may  coexist  within  the  bubble. 
The  variable  w has  therefore  been  introduced  to  represent  the  fraction  of 
the  hydrogen  which  has  actually  been  involved  in  the  reaction.  In  the  fuel 
lean  case  w will  be  -1.0,  while  in  the  fuel  rich  case  w will  be  - .42y  . 
The  factor  6 is  the  ratio  of  water  formed  to  hydrogen  Involved  and  has 
been  introduced  to  take  dissociation  into  account.  An  approximate  relation 
was  formulated  for  B(T,R)  as  a function  of  temperature  for  equivalence 
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Empirical  relation 

Totol  combustion  of 

0 homogeneous  mixture 


The  function  Is  normalized  to  the  mole  fraction  of 
1500*K  and  can  be  represented  as  follows: 


ratios  .5  ^ R ^ 4 . 
water  formed  at  T2 


where 


k 


1277.8 


826. 4|r-  1 
R + 2.375 


(12) 


(13) 


The  agreement  between  this  numerical  approximation  and  the  equilibrium 
calculations  of  Reference  9 Is  Illustrated  In  Figure  6.  Below  1500”K, 
virtually  all  of  the  hydrogen  Involved  In  the  reaction  will  be  converted  to 
water  and  0 will  equal  one. 

The  expansion  ratio  Is  related  to  the  simplified  chemical  change 
via  the  equation  of  state 


„ V2  ^2 

"2^STD  "^STD 

where  m Is  the  molecular  weight,  the  ambient  temperature,  and 

are  the  mole  fractions  derivable  from  Eq.  11. 

The  temperature  T2  of  the  combustion  products  Is  found  by  assuming 
that  the  combustion  and  expansion  process  Is  adiabatic  so  that 

“2(’^2)  " ‘*l(‘^?)  “ ‘*i('^STd)  ^ 

where  H2  la  the  enthalpy  of  the  combustion  products  and  Is  the 

enthalpy  of  the  bubble  before  Ignition.  The  molar  specific  heat  for  each 
species  can  be  expressed  as  a polynomial*** : 


(16) 


c 

Pi 


L,  + L,  4)  + L- 
ii  3^ 


+ L (()»  + L 4)-* 


where  4>  “ T/1000  and  the  coefficients  are  found  from  curve  fits  of  Indi- 
vidual species  In  the  JANNAF  tables** . Since 


where  k ■ can  be  expressed  as  a polynomial  as  well  and  Iterated 

with  respect  to  temperature  until  Eq.  IS  Is  satisfied  for  a given  w and  y. 
The  temperature  T2  can  then  be  Iterated  with  respect  to  w or  y to 
obtain  a solution  which  will  satisfy  both  Eq.  14  and  Eq.  15. 

For  the  case  of  complete  mixing,  hydrogen  and  oxygen  are  allowed  to  bum 
until  one  species  Is  depleted.  In  this  case  w becomes  unity  when  the 
bubble  Is  fuel  lean  and  w - .42y  when  fuel  rich,  where  y (v-Vg]/Vg  , 
the  volumetric  ratio  of  entrained  air  to  hydrogen.  This  simplistic  approach 
leads  to  an  expansion  time  history  which  differs  considerably  from  the 
available  data  as  can  be  seen  by  comparing  the  dashed  curve  with  the  Gorev 
hydrogen  data  In  Figure  5. 

There  are  several  reasons  why  the  data  do  not  follow  the  complete  mixing 
curves.  Since  sufficient  air  Is  entrained  for  complete  combustion  (stochlo- 
metrlc  mixture  ratio)  In  such  a short  time  (-.5  seconds  for  aim  diameter 
hydrogen  bubble).  It  Is  unlikely  that  the  hydrogen  and  air  are  mixed  on  a 
microscopic  scale.  In  the  fuel  lean  case  for  example.  It  Is  likely  that 
pockets  of  gas  exist  In  which  the  mixture  ratios  are  above  or  below  the  flam- 
mability limits.  Therefore,  not  all  of  the  hydrogen  will  be  burned.  Thus, 

In  the  fuel  lean  case  the  y of  Eq.  11  Is  determined  by  the  volume  of  en- 
trained air  and  the  value  of  w Is  decreased  according  to  the  Iteration 
scheme  discussed  above  until  the  appropriate  expansion  Is  obtained. 
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When  Che  mixture  Is  fuel  rich,  y cannot  be  found  directly  from  Che 
expanded  volume  because  Che  fuel  rich  mixture  apparently  reacts  with  Che 
surrounding  air  as  well  as  with  that  which  has  been  entrained*^.  In  this 
case,  y Is  Increased  beyond  the  amount  of  entrained  air  to  Include  the 
surrounding  air  Involved  In  the  combustion.  Assuming  chat  all  of  the  oxygen 
Is  depleted,  a value  of  y Is  found  which  satisfies  the  expansion  given  by 
Eq.  lA,  thus  determining  w6  as  well.  For  extremely  fuel  rich  mixtures 
%diere  the  expansion  after  Ignition  appears  to  be  very  large,  w Is  set  to  1 
and  y Is  Increased  beyond  Che  point  where  all  of  the  hydrogen  has  been 
depleted.  This  procedure  accounts  for  the  entrainment  of  air  which  does  not 
bum  with  the  hydrogen,  but  does  contribute  to  the  final  volume. 


Using  Che  above  analysis,  the  percentage  of  hydrogen  burned  and  the 
temperature  of  the  combusted  mixture  as  a function  of  time  after  release  are 
presented  In  Figures  7 and  8 for  bubbles  of  Initial  diameters  of  1 meter  and 
20  meters,  respectively.  The  percentage  of  hydrogen  In  the  bubble  prior  to 
Ignition  Is  also  shown  on  the  abscissa.  From  the  figures  It  can  be  seen 
that  the  percentage  of  hydrogen  burned  la  always  greatest  If  Ignited  Immedi- 
ately after  release.  The  highest  temperature  Is  achieved  when  the  bubble 
Is  Ignited  just  before  It  reaches  a stoichiometric  mixture  ratio. 


4.  DETERMINATION  OF  IR  RADIATION 


The  radiation  emitted  from  hydrogen  bubbles  after  combustion  is  also 
of  Interest,  both  as  a means  of  interpreting  the  readings  obtained  from 
various  sensors  and  as  testing  models  of  mixing  and  cooling  by  comparison 
with  existing  radiation  data.  Two  cases  of  interest  can  be  identified: 
(a)  the  radiation  immediately  after  combustion  as  a function  of  ignition 
tloie  and  (b)  the  radiation  time  history  of  the  bubble  after  ignition. 

An  isothermal  model  is  used  to  compute  the  radiation.  In  this  model 
the  Infrared  radiant  Intensity  can  be  expressed  as 


(17) 


in  watts/str/y  where  is  the  effective  emlsslvity  over  the  bandwidth  of 

interest,  ;^(T)  is  the  blackbody  function  for  wavelength  \ and  temper- 
ature T , and  > nr^  is  the  effective  area  of  the  radiating  bubble^*. 

The  emlsslvity  is  given  as: 


Ex  - 1 - TT  Ti 
^ 1-1 


(18) 


where  the  transmissivity  of  each  species  1 is: 


-Wi/d, 


(19) 


end  ■ ^1^1  weakline  approxlmetlon.  k^  is  the  absorption 

coefficient  in  the  2.8  micron  region  and  is  a function  of  temperature. 
is  the  optical  depth  and  is  expressed  as 


*1  “ ‘'lP»^eff^STD/’^2 


(20) 
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upper  radiation  limit.  VHien  the  mixture  ratio  of  the  homogenous  bubble  becomes 
stoichiometric,  however,  Its  radiation  signature  will  be  slightly  higher  than 
that  of  the  unmlxed  bubbles  since  dissociation  effects  are  less  severe.  Bubbles 
with  Initial  diameters  of  1,  5,  10,  and  20  meters  are  examined.  The  relative 
effects  of  Initial  upward  velocity  and  crosswind  on  the  Initial  radiation  from 
hydrogen  bubbles  of  1 and  20  meters  In  diameter  are  shown  In  Figures  10  and  11. 
The  effects  are  small  compared  with  the  radiation  dependence  on  Initial  bubble 
radius. 

After  Ignition  the  bubble  continues  to  rise.  As  it  mixes  with  cold  air 
the  radiation  drops  rapidly.  The  moles  of  air  entrained  at  any  time  t 
after  Ignition  can  be  found  using  the  Shul  and  Weyl  model  with  appropriate 
Initial  conditions.  The  temperature  of  the  bubble  Is  again  determined  by 
iteration  on  the  enthalpy  polynomials  from  JANNAF"  until  the  adiabatic 
expression  (Eq.  15)  Is  satisfied  for  the  given  ratio  of  entrained  air  to  hot 
combustion  products.  This  method  assumes  that  the  bubble  Is  thoroughly  mixed 
after  combustion.  The  volume  of  the  bubble  from  which  we  derive  the  effective 
path  length  Is  recalculated  from  the  equation  of  state  at  the  new  temperature, 
thus  taking  account  of  the  temperature  and  species  dependence  of  the  specific 
heat. 

An  additional  effect  that  must  be  accounted  for  In  the  calculation  of 
the  effective  temperature  Is  radiative  loss.  Some  of  the  energy  In  the  hot 
bubble  will  be  radiated  away  according  to  the  Stefan-Boltzmann  law,  J <■  eOT", 
where  c Is  the  engineering  emlsslvlty* The  temperature  drop  due  to 
radiative  loss  In  a time  At  Is: 

AT  - JAt/nCp  (23) 

where  the  temperature  drop  for  each  Increment  At  Is  small  with  respect  to 
the  bubble  temperature.  Recombination  of  the  reacting  species,  however,  would 
tend  to  maintain  the  bubble  near  the  flame  temperature  until  the  energy  of 
recombination  had  been  radiated  away. 
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Figure  10.  Effect  of  crosswind  and  upward  velocity  on  initial 

radiant  intensities  for  1 meter  diameter  H2  bubble. 
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Figure  11.  Effect  of  crosswind  and  upward  velocity  tTj  on  Initial 
radiant  Intensities  for  20  meter  diameter  H2  bubble. 


A. 2 Radiation  from  Combusted  Methane-Oxygen  Bubbles 


In  the  GEST  experiments  the  radiation  time  history  of  the  methane-oxygen 
bubble  after  Ignition  was  measured.  In  order  to  test  the  validity  of  the 
simplified  radiation  model  discussed  above,  radiation  estimates  for  the  hot 
combustion  products  were  obtained  by  applying  the  radiation  model  to  stoichio- 
metric methane-oxygen  bubbles  and  the  results  compared  to  the  GEST  data.  In 
our  model  the  atmosphere  Is  assumed  to  contain  .03%  CO2  by  volume  with  a 
humidity  of  20%  (1.34%  water  by  volume).  Figure  12  presents  theoretical 
curves  calculated  on  the  basis  of  three  different  entrainment  coefficients 
and  neglecting  the  effects  of  radiative  loss  and  recombination.  Again,  an 
entrainment  coefficient  of  .1-.12  Is  suggested  since  radiation  calculations 
on  the  basis  of  a larger  a fall  off  much  too  rapidly. 

As  a becomes  smaller  (a < .1),  the  form  drag  neglected  In  our  analysis 
may  become  Important.  If  form  drag  la  not  negligible,  the  entrainment 
coefficient  for  the  methane-oxygen  bubble  could  be  somewhat  smaller  (a “.OS) 
and  still  result  In  a radiation  prediction  similar  to  that  produced  by 
larger  a when  form  drag  effects  are  Included. 

Radiative  loss  and  recombination  effects  are  treated  In  Figure  13. 

Neither  effect  greatly  changes  the  radiation  time  history.  Although  a 10 
meter  diameter  stoichiometric  methane-oxygen  bubble  Immediately  after  Ignition 
radiates  at  the  rate  of  2.02  x 10^  watts/str/sec  corresponding  to  a tempera- 
ture loss  rate  of  lA78”K/sec,  the  cooler  bubble  Is  denser  at  each  point  In 
time  and  mixes  less  rapidly  with  the  air.  In  addition,  the  rate  of  temperature 
decrease  due  to  radiative  loss  drops  to  about  10”K/sec  In  .8  second.  Inclusion 
of  recombination  reduces  the  effect  of  radiative  loss  even  further. 

In  the  GEST  experiments,  spectra  of  the  radiating  bubble  In  the  region 
from  2400-6400  cm~^  were  recorded  at  .11  and  5.15  seconds  after  Ignition.  In 
Figure  12  the  data  at  .11  second  are  compared  to  the  SPECRA  band  model  dis- 
cussed above.  The  Input  parameters  were  bubble  radius,  mole  fraction  of 
radiating  species,  and  temperature.  Since  values  for  these  parameters  were 
taken  from  our  simple  model,  the  comparison  provides  an  additional  test  of 
the  model. 
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Radiation  time  history  of  a stoichiometric  CH^-02  bubble  after  ignition,  X - 2.8 
Dependence  on  entrainment  coefficient. 
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Figure  13.  Radiation  time  history  of  a stoichiometric  CH^-02  bubble  after  ignition,  X ■ 2,8 
Effect  of  radiative  loss  and  recombinstlon. 


Figure  14.  Spectral  radiant  Intensity  for  CHa-Oo  bubble  at  .11  second  after  Ignition 


The  data  have  been  corrected  for  atmospheric  absorption,  but  In  some  wave- 
length regions  the  spectral  radiant  intensity  was  so  low  due  to  the  severity 
of  the  attenuation  that  accurate  adjustments  could  not  be  made.  The  SPECRA 
calculation  is  somewhat  higher  than  the  data,  but  well  within  the  scatter  of 
the  data  (see  Figures  12  and  13).  The  same  variation  with  X is  evident  in 
both  data  and  theory  except  in  those  wavelength  regions  where  the  data  were 
severely  attenuated. 

4.3  Radiation  Predictions  for  Hydrogen  Bubbles  After  Ignition 

In  Figures  15-18,  the  model  neglecting  radiative  loss  and  recombination 
is  applied  to  hydrogen  bubbles  1 and  20  meters  in  diameter.  Figures  15  and  17 
show  altitude-time  histories  for  different  ignition  times.  The  dashed  line  is 
the  rise  time  history  of  an  uncombusted  hydrogen  bubble.  As  the  uncombusted 
hydrogen  bubble  rises  and  mixes  with  the  air,  it  rapidly  drops  below  the 
flanmabllity  limit.  By  the  time  a hydrogen  bubble  20  meters  in  diameter  has 
risen  to  a height  of  40  meters  (10  seconds),  it  cannot  be  Ignited.  A 1 meter 
diameter  bubble  reaches  this  limit  in  only  2.2  seconds  at  a height  of  2 meters. 
The  rapidity  of  this  process  is  partially  due  to  the  unmlxedness  of  the  bubble 
before  ignition. 
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Figure  15.  Altitude  time  history  at  various  ignition  points  fo 
1 meter  diameter  hydrogen  bubble  starting  from  rest 


Altitude  tlae  history  at  various  Ignition  points  for 
20  aeter  disaster  hydrogen  bubble  starting  froa  rest 
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I Figurt  18.  Radiation  tlM  history  st  various  Ignition  points  for  a t 

20  Mtar  dlaactar  hydrogen  bubble  starting  froa  rest.  I 
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5.  SUMMARY 

The  motion  of  buoyant  gas  clouda  has  been  examined  using  the  simplified 
model  found  In  References  1 and  5.  The  effect  of  a crossvlnd  on  bubble 
motion  la  also  treated  but  Is  found  to  have  negligible  effect.  Comparison 
of  the  model  with  data  on  small  hydrogen  bubbles  at  ambient  temperature, 
hot  combusted  CH^-02  clouds,  and  thermals  caused  by  nuclear  explosions  Indi- 
cates that  the  model  Is  applicable  to  a wide  variety  of  bubble  types  with 
the  appropriate  choice  of  entrainment  coefficient.  Since  there  Is  not 
sufficient  data  at  this  time  to  concluslvo.ly  determine  the  entrainment 
coefficient  suitable  for  large  H2  bubbles,  the  entrainment  coefficient 
measured  by  Gorev  et  al.  for  small  hydrogen  bubbles  (a  ■ .41)  Is  used. 

Gorev  et  al.  Ignited  H2  bubbles  at  various  altitudes  and  measured  the 
volume  of  the  combustion  products  at  these  different  mixture  ratios.  The 
degree  of  expansion  predicted  by  total  combustion  of  a homogenous  mixture 
differs  considerably  from  the  data  since  the  bubble  Is  never  mixed  on  a 
microscopic  scale  and  bums  with  the  surrounding  air  when  fuel  rich.  A 
seml-emplrlcal  theory  was  developed  whereby  parameters  of  the  combustion 
process  such  as  percentage  of  H2  burned  and  temperature  after  Ignition  could 
be  obtained  from  a model  which  maintained  the  observed  expansion.  It  was 
found  that  the  temperature  of  the  cloud  after  Ignition  Is  highest  If  the 
bubble  Is  Ignited  just  before  It  reaches  a stochlometrlc  mixture  ratio.  In 
addition,  the  more  fuel  rich  the  mixture,  the  greater  the  percentage  of 
hydrogen  which  Is  converted  Into  water. 

A lower  and  upper  limit  to  the  radiation  Is  obtained  by  applying  the 
radiation  model  developed  by  Flshbume**  to  bubbles  In  which  the  combustion 
products  and  unreacted  species  are  either  thoroughly  mixed  or  unmlxed, 
respectively.  The  time  history  of  the  radiation  after  Ignition  as  the  bubble 
cools  Is  also  examined  and  compared  with  data  on  combusted  CH^-02  clouds. 
Radiative  loss  and  recombination  were  also  Included,  but  seem  to  have  only 
a small  effect.  The  entrainment  coefficient  Is  a much  more  sensitive 
parsmeter.  The  model  was  spplled  to  H2  bubbles  after  Ignition  as  «wll.  A 
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20  meter  diameter  bubble  was  seen  to  have  a peak  radiation  of  10^  watts/str/y 
If  Ignited  within  several  seconds  after  release  while  still  In  the  fuel  rich 
domain.  However,  after  10  seconds  and  at  an  altitude  of  40  meters,  the  bubble 
becomes  so  fuel  lean  that  It  can  no  longer  be  Ignited.  The  smallest  bubble 
considered  (a  diameter  of  1 meter)  cannot  be  Ignited  after  only  2.2  seconds 
or  above  an  altitude  of  2 meters.  It  thus  reaches  Its  peak  radiation  of 
10^  watts/str/y  only  If  Ignited  within  .6  second  after  release.  The  radiation 
of  the  hot  combustion  products  drops  off  very  rapidly  after  Ignition.  Even 
for  the  20  meter  diameter  bubble,  the  maximum  time  In  which  a significant 
radiation  signature  Is  sustained  Is  8 seconds. 
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